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ABSTRACT 

We describe a new tool for studying the structure and physical characteristics of ultracompact 
AGN jets and their surroundings with /ias precision. This tool is based on the frequency dependence 
of the light curves observed for intra-day variable radio sources, where the variability is caused by 
interstellar scintillation. We apply this method to PKS 1257-326 to resolve the core-shift as a function 
of frequency on scales well below ~ 12 /ias. We find that the frequency dependence of the position of 
the scintillating component is r ex ^- 01 ± - 24 (99% confidence interval) and the frequency dependence 
of the size of the scintillating component is d cx ^-°- 64 ± 006 . Together, these results imply that the 
jet opening angle increases with distance along the jet: d cx r nd with > 1.8. We show that the 
flaring of the jet, and flat frequency dependence of the core position is broadly consistent with a model 
in which the jet is hydrostatically confined and traversing a steep pressure gradient in the confining 
medium with p cx r~ Up and n p > 7. Such steep pressure gradients have previously been suggested 
based on VLBI studies of the frequency dependent core shifts in AGN. 

Subject headings: galaxies: jets — techniques: high angular resolution — quasars: individual 
(PKS 1257 - 326) — scattering 



1. INTRODUCTION 

The brightest, most compact feature of an AGN 
jet, the "core", is identified with the part of the 
j et at which the optical depth (r u ) is of order unity 
(Blandfo rd fc Konigll Il979h . and is often referred to as 
the t v = 1 surface, or photosphere. Due to positional 
variation of the opacity in the jet and/or surround- 
ing medium, the position of the t v = 1 surface is fre- 
quency dependent, and th erefore, so to o is the absolute 
position of the cor e (eg. iKo nigl 1981; Lobanovl 119981 : 
iKovalev et all [20081 : ISokolovskv eTall [20TTh . The fre- 
quency dependent position of the core is referred to sim- 
ply as the core shift. 

The core shift effect provides an observational tool 
with which to investigate the structure and physical con- 
ditions in parsec-scale AGN jets. Moreover, modelling 
the effect may provide information about the confine- 
ment mechanism and pressure gradients in the external 
medium. The core shift effect is also relevant to the quest 
for high precision absolute astrometry for the Interna- 
tional Celestial Reference Frame, as it can introduce a 
significant offset in positions det ermined using group de- 
lay measurements (jPorcasll2009[ ). 

The magnitude of the core shift between 2.3 and 
8.4 GHz is typically of order a f ew hundred /ias 
or le ss (10'Sulliva n & Ga buzdal 120091 : ISokolovskv et al.l 
2011; Pushkarev et al. 2012), and therefore detecting this 
effect requires very high accuracy registration of im- 
ages at two or more frequencies. Despite the techni- 
cal challenges, the core-shift effect obtained from VLBI 
imaging has been re ported for an ever- increasing num- 
ber o f radio galaxies ( Marcaide & Shapiro 1984; Lob anovl 
1998; Kovalev et al. 2008; O'Sullivan & Gabuzda 2009; 



1 ARC Centre of Excellence for All-Sky Astrophysics (CAAS- 
TRO) 

2 Now at Max-Planck-Institute fur Radioastronomie, Auf dem 
Hiigel 69, 53121, Bonn, Germany 



ISokolovskv et al.ll20lH IF ushkarev et al.ll2012f ). More re- 
cently, |Kudryavtseva et al] (|2011l ) have employed an in- 
direct method to measure the core shift effect based on 
frequency-dependent time lags of flares observed using 
single-dish data spanning several years. 

The frequency dependence of core position is typically 
assumed to follow a power-law of the form r cx v~ x l kr . 
In many sources for which core shifts can be measured 
with VLBI imaging, the absolute core position varies 
approxi mately with the inverse of the frequency (i. e. 
k r = 1) (jO'Sullivan fc Gabuzdal 12001 ISokolovskv et al.l 
|2011[ ). This situation is consistent with the standard 
model of a conical jet in which the plasma is in a state 
of equipa rtition between particle an d magnetic energy 
densities (Blandf ord fc Konigll 1 19791 ). However, values 
of k r much greater than unity are observed in some 
sources, which may be due to free- free absorption in the 
immediate vicinity of the jet, or due to rapid changes 
in pressure in the exte rnal medium i f hydrostat i c con - 
finement is important (Lobanov 1998). Lobanov (1998) 
has shown that while k r ~ 1 at large distances down- 
stream from the black hole, the value of k r inc reases 
towards the jet base. iKudrvavtseva et al.l (|2011f ) have 
shown that the value of k r is time-dependent, and corre- 
lated with flux density. Finally, the pc-scale jet of M87 
is observ ed to deviate from a conic al geometry near to 
the core (As ada fc Nakamural 120121 ). Further investiga- 
tion into the frequency dependence of core position is 
therefore warranted, and highly relevant to the study of 
ultracompact AGN jets. 

Here we present a potentially powerful new method 
for the study of ultracompact jets in AGN, which enables 
simultaneous measurement of the core shift effect and jet 
geometry to very high precision. This new technique uses 
auto- and cross-correlation analysis of multi-frequency 
light curves of a rapidly scintillating AGN to measure 
the frequency dependence of the position and size of the 
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scintillating component. 

In section [2] we present the observations and data anal- 
ysis. In Section [3] we discuss the mathematical formula- 
tion of the auto- and cross-correlation analysis, and de- 
rive the frequency dependent source position and size for 
PKS 1257-326. In Section [4] we discuss the implications 
of our findings, and model the jet in terms of a hydrostat- 
ically confined jet traversing a steep pressure gradient. 
Finally, in Section [5] we present our conclusions. 

2. OBSERVATIONS AND DATA CALIBRATION 

PKS 1257-326 was observed at the ATCA for ten 
hours on 15 January 2011, with two 2 GHz bands, cen- 
tred on frequencies of 5.5 and 9.0 GHz. The output data 
included all four polarisation products and 2048 spectral 
channels each 1 MHz wide in each of the two bands. Flag- 
ging and calibration of the data were performed using the 
Miriad software package. The ATCA primary calibrator 
PKS 1934—638 was used to correct the overall flux den- 
sity scale and the spectral slope. In order to solve for the 
bandpass and to correct gain amplitudes as a function of 
time and pointing for each antenna, the secondary cali- 
brator PKS 1255-316, only 1° from PKS 1257-326, was 
observed for 1 minute approximately every 20 minutes, 
interleaved with observations of the target source. Phase 
self-calibration assuming a point source model was per- 
formed with a short (10 s) solution interval. After initial 
calibration the data were split into 128 MHz sub-bands 
for further analysis. 

At least 98% of the total flux density of PKS 1257-326 
is unresolved with the ATCA, and there is no significant 
confusion in the field at frequencies above 4 GHz. There- 
fore, to obtain the light curves of the source variations we 
averaged the real part of the calibrated visibilities over 
all baselines and frequency points within each 128 MHz 
band. The relative stability of the calibration as a func- 
tion of time and frequency within each 2 GHz band is esti- 
mated to be ~ 1% or better, based on the PKS 1255—316 
data. The bandpass is observed to be stable over the 
duration of our observations in the 5 GHz band, but 
there are small frequency- and elevation-dependent gain 
variations across the higher frequency band, which were 
corrected with a time-dependent bandpass solution de- 
rived from PKS 1255— 316. Although the frequency de- 
pendence of the primary calibrator PKS 1934—638 is 
well known, archival data on the secondary calibrator 
PKS 1255—316 shows it to be variable by up to ~ 50% 
on timescales of months to years. Hence there is a 
small uncertainty in the spectral slope correction for the 
9 GHz band, due to the variations with time and point- 
ing, and the fact that the spectrum of PKS 1255-316 is 
not known a priori. In any case, the average spectrum 
of PKS 1257—326 is relatively smooth across the entire 
range of frequencies, suggesting that the calibration is 
accurate. Moreover any residual constant offsets in the 
flux density scale which may be present between different 
frequencies have no effect on the cross-correlation analy- 
sis presented in this paper. 

Figure [T] shows the large, rapid intra- hour variations 
exhibited by PKS 1257-326. 

3. DERIVATION OF JET STRUCTURE 
3.1. Interpretation 



The rapid fluctuations observed in the centimetre 
wavelength flux density of the intra-hour variable quasar 
PKS 1257—326 are due to interstellar scintillation (Big- 
nail et al. 2003, 2006). This is established from the mea- 
surement of a time delay in the arrival time of the vari- 
ations between telescopes separated by several thousand 
kilometers. The timescale of the variations also under- 
goes an annual modulation due to relative motion of the 
Earth about the Sun, which in turn moves relative to the 
interstellar medium responsible for the variations. 

Inspection of Figure [T] reveals that there is a time offset 
in the arrival time of the intensity variations between dif- 
ferent frequencies, with the variations at high frequency 
leading those at lower frequencies. This behaviour is con- 
sistent with observations in previous epochs. Bignall et 
al. (2003) reported that variations at 4.8 and 8.6 GHz are 
closely correlated, and there is a systematic time delay 
between the variations at these two frequencies. The 
magnitude of the time delay was observed to follow an 
annual cycle which is not identical to the annual cycle in 
variability timescale but, it was argued, can be explained 
on the basis of it, as we discuss below. 

PKS 1257-326 was monitored with the ATCA at 4.8 
and 8.6 GHz at 19 epochs between 2001 February and 
2002 April. Typical observations were over a 12 hour 
period, and 6 epochs covered 2 x 12 hours in 48 hour 
sessions as part of a multi-source monitoring program. 
The minimum duration of each light curve is ~ 10 times 
the length of the characteristic timescale. In every one 
of these epochs, the time delay between 4.8 and 8.6 GHz 
has the same sign, with 8.6 GHz variations always lead- 
ing. Moreover, a clear annual cycle is observed in the 
two-frequency time delay, with the longest delays being 
observed from late July to mid- August. Such an annual 
cycle is expected for a core-shift which remains stable 
over the course of the year, and the observed annual cy- 
cle is well modelled by such a shift on a scale of order 
10/ias (Bignall et al. 2003), although the precise magni- 
tude and direction of the core shift could not be uniquely 
determined. These data provide strong evidence that the 
core shift effect dominates over any refractive effects or 
"jitter" in the ISS pattern. 

We argue that the temporal offset in the present data 
arises as a direct consequence of an angular offset be- 
tween two compact components within the scintillating 
source. The effect may be understood as follows. When 
an angular separation, 0, is present between two com- 
ponents, this results in a spatial displacement of their 
respective scintillation patterns across the plane of an 
observer by an amount DO, whe re D is distance be- 
tween the observer and the scattering material (Little & 
Hewish 1966). Since the scintillation patterns are in mo- 
tion across this plane with some velocity v, the result 
is a separation in the arrival time of the scintillations 
associated with each component. In the present case, a 
displacement in the lightcurves between closely-spaced 
frequencies arises because there is an angular offset in 
the image centroids between the respective frequencies. 
For any pair of frequencies, the time delay is, in terms of 
the centroid offset (see Appendix |A|) , 
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Fig. 1. — PKS 1257— 326 flux density measurements as a function of time, plotted with 1-minute averaging, showing each 128 MHz band 
in a different colour. Frequency increases from top to bottom. The heavy black points show the average of each 2 GHz wide-band dataset, 
corresponding to the centre frequencies of 5.5 and 9 GHz. 



where R is the anisotropy ratio of the scintillation pattern 
and S = (cos /?, sin f3) is the direction of its major axis, 
which we measure with respect to the RA axis. The scin- 
tillation parameters have been derived from annual cycle 
and two-station time-delay measurements, and are given 
in Table 1. It is evident that this delay is modulated both 
by the annual cycle in the magnitude of v and by changes 
in the angle of the velocity vector with respect to 0; this 
latter effect causes the annual cycle experienced by At 
to differ from the annual cycle in scintillation velocity. 

Phase gradients in the ISM may, in principle, also cause 
temporal offsets of lightcurves as a function of frequency 
in a scintillating source. However, the offset observed 
here is difficult to attribute to such an extrinsic cause 
for several reasons: (i) the sense and magnitude of the 
delay is constant throughout the dataset; upon dividing 
the dataset in two halves (in time) and deriving time 
offsets based on these two halves separately, we find 
the same offsets to within the margin of error of the 
estimates, (ii) The delay is observed over a timescale 
of 10 hours, whereas refractive phase gradients in the 
ISM in the regime of weak scintillation for a Kolmogorov 
spectrum of phase inhomogeneities would occur on the 
timescale associated with the scintillations, and the time 
offset should converge to zero as the average is performed 
over an increasing number of scintles. Any small jitter in 
the offset between individual scintles appears to be dom- 
inated by the systematic offset, (iii) An annual cycle in 
the time offset is reported by Bignall et al. (2003), in- 



dicating that the offset persists on a timescale of greater 
than a year. 

3.2. Time delay measurement 

To determine the relative time delay between each pair 
of lightcurves, I(t,v\) and /(t,^), we computed the 
cross-correlation function, 



C(At;i/i,i/ 2 ) 



<[/(*', i/i) - I(t, z/i)] [lit' + At, v 2 ) - lit, v 2 )}) 



y/vax[I(t, ^i)]var[/(t, v 2 ) 



(2) 



A peak in the cross correlation at positive delay, At, 
indicates that the fluctuations at frequency v 2 precede 
those at v\ . We fitted a gaussian of the form, 



C{t) = Aexp 



(t-to) 2 



B 2 



(3) 



to the inner part of the cross-correlation function, 
C(At; v\, v\) (equivalent to the auto-correlation func- 
tion), to obtain an estimate of the time delay between 
each frequency- lightcurve pair and its associated error. 
An example cross-correlation function and its associated 
fit is shown in Fig. [2] Typical errors in the estimated de- 
lay are 50 s. The derived delays as a function of v\ and v 2 
are shown in Fig. 03 The estimated 50 s uncertainties in 
the delay between each frequency pair are derived from 
least-squares fitting of a gaussian to the peak of the de- 
lay, thus the errors are directly related to the width of 
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TABLE 1 

The parameters of the scintillations in PKS 1257-326. 



Parameter^ 


Symbol 


B06 value 


W09 value 


screen distance 


D 


10 pc 


10 pc 


scintillation velocity, v, 


Vol 


54.12 kms" 1 


SS.Tkms- 1 


on observation date 


V S 


-0.05 kms" 1 


Onions -1 


anisotropy ratio 


R 


12 


6 


anisotropy orientiation 


(3 


-35° 


-33° 



a Parameters are derived on the basis of the annual cycle in the vari- 
ability timescale (Bignall et al. 2003) and measurements of the time 
delay observed between two stations (Bignall et al. 2006, henceforth 
B06). An alternate fit to the scintillation data provided by Walker, 
de Bruyn & Bignall (2009) (henceforth W09) is also listed. The scin- 
tillation velocity quoted here is the addition of Earth's relative to the 
Sun velocity on the observation date and the velocity of the scattering 
medium relative to the Sun. 



the peak in the cross-correlation function^. These un- 
certainties are in turn derived from the least squares fit 
to the cross-correlation function; the errors in the indi- 
vidual points in the cross-correlation function are dom- 
inated, at low time lags, by Poisson errors associated 
with the number of independent measurements of the 
cross-correlation measureable from the lightcurve. This 
error is indicated by the scatter between points in the 
CCF. Fitting to the CCF yields a typical formal error 
of ~ 20 s. However the high degree of cross correlation 
between the lightcurves in our observations means that 
the error in the time delay estimate for any given fre- 
quency pair is not completely independent from the time 
delay estimate of any other frequency pair. It is neces- 
sary to take the cross-correlation into account because it 
becomes an important factor when assessing the signifi- 
cance of the fit to the frequency dependence of the time 
delay using a least-squares approach (see below), which 
is the primary reason to estimate the error in the time 
lag. Failure to take into account this interdependence in 
time lag estimates would lead to an overestimate of the 
significance of the fit to the frequency dependence of the 
time shift. It was found empirically, from examination 
of the reduced chi-squared in the fit procedure that this 
cross-correlation is taken into account in the error anal- 
ysis with errors that are 2.5 times larger than the formal 
error. 

In most models of jet structure the angular offset be- 
tween the jet base and the centroid of the emission at 
a frequency v is fitted to a power law, AO = Av~^ 
where A and ( are constants to be determined. Since 
the time delay is linearly proportional to |0|, we fitted a 

function of the form At (^1,^2) = — 1^2^) to the 

delays. The best-fitting parameters are ( = 0.10 and 
K = (1.6 ± 0.5) x 10 4 sMHz^, with the 99% confidence 
interval of ( extending over the range [—0.14,0.34]; the 
best-fitting surface is plotted over the data in Fig. [30. 

3 This error is significantly smaller than the time delay errors 
that might be estimated on the basis of the data presented in Big- 
nall et al. (2003). As is evident from Fig. 2 in Bignall et al. (2003), 
the data in most of those epochs contain considerably few scintles 
from which to estimate the time lag. Thus one would expect the 
errors to be larger relative to the 2011 data. The dataset obtained 
on 4 Jan 2001 contains a comparable number of scintles, but here 
the temporal sampling of the lightcurve was too sparse to estimate 
a time delay. 

4 The 99% confidence limits are calculated, in the conventional 
way, from the change in the sum of the squares of the residuals 



The best-fitting exponent indicates a core-shift depen- 
dence on frequency that is much shallower than the value 
C = 1 typically found in other quasars on the basis of 
VLBI measurements (e.g. 10 'Sullivan fc Gabu zda 2009; 
iSokolovskv eTaII[20TTh . 

We also measure the timescale of the scintillations, 
which is derived from the parameter B in eq. (j3j) (and 
with v\ = v^). Figured] shows that the timescale is well- 
fit by a relation scaling as ^-°- 642 ± 006 - a fit using a 
broken power law reveals the both the 4-5 and 8-9 GHz 
timescales follow the same scaling with frequency within 
the (small) margin of error. 

In the regime in which the angular size of the source, 
# src , exceeds the angular size of the Fresnel scale at the 
distance of the scattering screen, Of = ^JcI^/kvD , the 
scintillation timescale is linearly proportional to the size 
of the scintillating component. In the opposite regime, 
#src < #f, the scintillation timescale is determined by 
DOf/v = rp/v, which scales as u~ - 5 . 

We wish to determine whether the timescale 
of the scintillations measures Of or src . 
The expected Fresnel crossing timescale is 
1.0 x 10 3 (D/10pc) 1 / 2 (^/5GHz)- 1 /2( v /54kms" 1 )- 1 s. 
The observed scintillation timescale at 5.0 GHz is 
1.4 x 10 3 s, so if the source is unresolved there must 
be an error in the nominal scintillation parameters; if 
the scintillation velocity is held at its nominal value 
one must have D > 19 pc, or if the screen distance is 
held at its nominal value one must have v < 40 kms -1 . 
The latter option is viable for the scintillation velocity 
derived by Walker, de Bruyn & Bignall (2009), so we 
conclude that it is plausible that the source is unresolved 
by the scintillations. However, the fact that the scaling 
of scintillation timescale is significantly different from 
v~ - 5 suggests that the source is at least partially 
resolved in any case. The smooth trend evident in both 
timescale and scintillating amplitude with frequency 
suggests that the source remains resolved over the entire 
frequency range 5-10 GHz. 

Although tangential to the objectives of this paper, we 
can, in addition, estimate the amplitude of the scintil- 



of the fit. Now, in practice, since the lightcurves are correlated 
between frequencies, the time lag estimate between adjacent fre- 
quency pairs are strictly not independent, and this influences the 
chi-squared estimate. However, this effect is treated by accounting 
for the cross-correlation in the estimates of the time delays for each 
frequency pair. 
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Fig. 2. — The cross correlation between the lightcurves at 8040 
and 4540 MHz, and the associated best-fit gaussian. The peak of 
the gaussian represents the location of the time-delay. The positive 
offset of the peak indicates that the variations at 8040 MHz lead 
those at lower frequency. 

lating component of the source. For an extended source 
of size SYC in the regime of weak scintillation caused by 
a Kolmogorov spectrum of turbulent fluctuations, the 
observed rms, (SI 2 ) 1 / 2 can be expressed in the form 
(Narayan 1992), 



(SI 2 ) 1 / 2 




(4) 



where (J) is the flux density of the scintillating compo- 
nent of the source and oc is the diffractive scale, 
which is determined by the properties of the interstellar 
turbulence (see Narayan 1992). In the regime of interme- 
diate scattering one has r^ « rp, with equality at the 
transition frequency, which likely occurs in the range 3- 
7 GHz based upon the modelling of Walker (1998). One 
then solves for (I) using the measurements of (SI 2 ) 1 / 2 . 
Assuming that the source size exceeds the Fresnel angu- 
lar scale, one determines # S rc/#F from the ratio £ SC int/£F- 
We performed a fit using the measured values of (SI 2 ) 
and tscint? and subject to the assumption that = 
at 4 GHz, to derive a rough estimate of the flux density 
of the scintillating component: 



/ v \0.5 



(5) 



We caution that the spectral index of the component de- 
rived here is only approximate because the expression for 
the modulation index in eq.(|4]) is only approximate at fre- 
quencies where r^if? ~ rp. A more sophisticated estimate 
would employ a more complicated approximation to the 
modulation index near the transition frequency and take 
into account the anisotropy of the scintillations. 

3.3. Overall source scale 

Although the frequency scaling of the core position and 
size are the primary results of our analysis, it is possible 
to roughly relate these measurements back to physical 
scales within the source. 

The Schwarzschild radius of the black hole at the cen- 
tre of PKS 1257-326 is R s = 2.95 x 10 12 M 9 m, where 
Mg = 10 9 M is the BH mass, which is estimated from 
measurements of the broad line emission width to be 
10 88 M© (D'Elia, Padovani & Landt, 2003). At the 
redshift of PKS 1257-326, z = 1.256, 1 /ias subtends 
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Fig. 3. — The time delay associated with each pair of lightcurves 
measured at frequencies v\ and h>2- The red surface corresponds 
to the best-fitting model of the form At = A(z/^ — is^**) ( see text 



for details). 




Frequency (GHz) 

Fig. 4. — The timescale of the scintillations as a function of fre- 
quency. The line represents the best-fitting model of the form 
r = Kv~^ (see text for details). Note that the intrinsic scat- 
ter between adjacent points is smaller than the error bars: this is 
because the lightcurves are highly correlated leading to a smaller 
scatter than is reflective of the errors. 

8.5 x 10" 3 pc, equivalent to 89/M 9 Schwarzschild radi0. 

The time delay is related to angular structure in the 
source using eq. (pQ) and the scintillation parameters in 
Table 1. This also depends on the angle, £, that 6 makes 
with the right ascension axis, which is unknown. For a 
time offset At, the expected amplitude of the angular 
offset is given by 



_ f (27.7 cos f ■ 
\ (43.9 cos f ■ 



-38.8sinf)" 
-61.8sin£)- 



x /ias, B06, 
Vas,W09, 



(6) 



where the two solutions denote the scintillation parame- 
ters found by Bignall et al. (2006) (B06) and Walker et 
al. (2009) (W09). For instance, the angular separation 
implied by the 520 s delay observed between the lowest 
(4540 MHz) and highest (9960 MHz) frequency bands of 
our observations is 19 /ias (12 /ias) for the scintillation 
parameters of B06 (W09) if £ = (i.e. if the angular sep- 
aration is aligned parallel to the right ascension axis). 



5 We use H 
throughout. 



TOkms- 1 Mpc-\ Q m = 0.27 and Q A = 0.73 
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This translates to a physical scale of 0.16 pc (0.10 pc) at 
the source. 

However, these estimates are subject to uncertainties 
in both the screen distance and the orientation of the 
separation £. For instance, doubling D above its nominal 
value of 10 pc would result in estimates of |0| lower by 
a factor of two. They thus serve only as a guide to the 
order of magnitude of scales which are probed by these 
measurements. 

In the same vein, we also estimate the angular sc ale of 
the source from the variability timescale, using eq. ([A 10)) 
in Appendix [A] to find the characteristic angular scale of 
the scintillation pattern, 



(9.96 GHz) 



-0.64 



D 
lOpc 





"506" 


1 /ias 


W09 { 



(7.) 



The minor axis of the scintillation pattern has a size 
STC /yfR while the major axis has a characteristic size 
of src VR. 

4. DISCUSSION 
4.1. Observational Constraints 

The foregoing analysis indicates that the centroid of 
the scintillating component has a frequency dependence 
of the form: r core oc ^-O-iiO- 24 . j n the terminology of 
iKonigll (|1981f ) and lLobanovl (|1998[ ). r core oc v~ x l k ^ we 
have k r > 3. We proceed under the assumption that the 
source is resolved, as suggested by the frequency depen- 
dence of the scintillation time-scale (Sect ion [32]). In that 
case, the jet diameter has a frequency dependence of the 
form: d oc ^- - 64 ± 006 . Taken together, these results im- 
ply d oc r™o re with rid > 1-8. The scenario implied by our 
analysis is illustrated in Figure [5] This is in contrast to 
the frequency dependence of the core position in VLBI 
studies: typically r core oc consistent with expecta- 
tions for a conical jet in which the part icle and magnetic 
energ y densities are in equipartition (So kolovsky et al.l 

I2QDJ- 
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Fig. 5. — A cartoon illustrating the constraints implied by our 
analysis. 



4.2. Interpretation in terms of hydrostatic jet 
confinement in a steep pressure gradient 

lLobanovl (|1998[ ) suggested that the changes in jet opac- 
ity observed for a sample of sources appear consistent 
with a self-absorbed jet propagating through a region 
with a steep pressure gradient of the form p oc r -6 . 
Here we consider such a model to explain the inferred 



frequency dependence of the core size and position in 
PKS 1257-326. Specifically, we consider whether such 
a model can simultaneously account for the increase in 
opening angle along the jet (d oc r™* re with rid > 1-8), and 
the flat frequency dependence of core position (r core oc 
v -i/k r ^ ^ifo fc r > 3) The relevant model equations, for 
a pressure profile p oc r~ Up , are (see Appendix |B|) 



n d = -r 



(n p /4)[(1.5 + a) (2 + m B ) - 1] - (1.5 + a)n 6 



2.5 + a 



5 I \ f m B 



■ ns 



Here vtib is a parameter determined by the magnetic field 
geometry as discussed in Appendix B, and a is the op- 
tically thin spectral index which directly relates to the 
electron energy distribution N(j) oc 7~( 2q! + 1 ), as op- 
posed to a core , which is the optically thick spectral index 
of the self-absorbed core. We assume a = 0.5, ra# = 1 
(i. e. a predominantly toroidal magnetic field configu- 
ration) and take n p = 7.2 so that rid = 1.8. In that 

case, the jet Lorentz factor goes as T/Tq = (r/ro) 1 ' 8 , 
so that for reasonable values of the jet viewing angle 
(5° < <j> < 40°) and initial Lorentz factor (2 < T < 10), 
we have —0.5 < n§ < 0.5, which implies: 



2.7 < k r 

0.7 < a CO re 



i$ 3.3, 
< 1.2. 



(8) 



There is excellent agreement between the predicted 
and observed frequency dependence of core position and 
radial dependence of core size. The agreement be- 
tween the predicted and observed core spectral index 
is poorer, but this is perhaps unsurprising since, as re- 
marked above, observational determination of the core 
spectral index is marred by uncertainties in the model 
used to derive it. We therefore suggest that this simple 
model of hydrostatic confinement is broadly consistent 
with the data, provided n p > 7.2. The required pressure 
gradient in the external medium (n p > 7.2) appears very 
steep, however, such pressure gradients have p reviously 
been suggested in studies of core shifts in AGN ([Lobanovl 
fl998h . 

Finally, we note that the frequency dependence of core 
position may also be influenced by free-free a bsorption 
in th e immediate environment of the jet (e. g. lLobanovl 
1998). However, we find no evidence for a significant 
rotation measure, and such an interpretation cannot ac- 
count for the flaring jet geometry, so that an additional 
mechanism, such as entrainment, would also be required 
to explain these results. 

5. CONCLUSIONS 

We have developed a new technique to probe the struc- 
ture of AGN jets on micro-arcsecond scales, by using 
interstellar scintillation to simultaneously determine the 
shift in the position of the AGN core as a function of fre- 
quency, and the frequency scaling of the core size, to high 
precision. This approach is amenable to sources which 
harbour compact (< 100 /ias) features. 

This method was applied to broadband observations 
of PKS 1257-326 with the Australia Telescope Compact 
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Array. The scaling of the core-shift is found to be re- 
markably shallow with frequency; the best fit to the posi- 
tion of the scintillating component in the source scales as 
r oc ^~ 010 , with the 99% confidence interval of the index 
extending over the range [—0.34,0.14]. This constrasts 
with previous VLBI studies which typically - though not 
exclusively - find r oc v~ x . The scaling of the jet size is 
also determined, based on the scaling of the scintillation 
timescale with frequency. This shows that the jet size 
scales formally as doc zy -o.64±o.oo6 j^. - g possible to de- 
termine the scaling of the core shift and jet diameter to 
high precision because they do not depend critically on 
complete knowledge of the properties of the scattering 
medium responsible for the scintillations. Determina- 
tion of the absolute physical scale of the core shift does, 
however, require knowledge of the scintillation param- 
eters, and we can only determine these quantities ap- 
proximately. The observed 520 s time offset between the 
scintillations at the the lowest (4540 MHz) and highest 
(9960 MHz) frequency bands of our observations implies 
an angular separation of > 12 /ias, for an assumed scat- 
tering screen distance of 10 pc. This translates to a phys- 
ical scale of ~ 0.10 pc at the source. This physical scale 
probed is an order of magnitude smaller than typical core 
shifts obtained with VLBI measurements. We further 
note that this technique easily detects the core-shift be- 
tween frequency pairs separated by only ~ 300 MHz, thus 
providing sub-microarcsecond resolution of the jet struc- 
ture. 

A major conclusion arising from our analysis is that 
that the often assume d frequency scaling of core posi- 
tion, r cor e oc v ~ x (e. g. iPushkarev et al.ll2012l ). may not 
be applicab l e to a ll so urces, in line with similar fi ndings 
of lLobanovi (|1998f ) and iKudrvavtseva et al.l (|2011l ). Fur- 
ther, our results hint at a physical difference between per- 
sistent intra-day variable (IDV) sources, and the broader 
population of AGN. 

To place these results in a physical context, we have 
explored a simple model based on a hydrostatically con- 
fined jet traversing a pressure gradient. The pressure 
profile implied by this model is steep: p cx r _Up ; n p > 7. 
Such steep pressure gradients have previously been sug- 
gested in VLB I studies of the frequency dependent core 
shifts in AGN (|Lobanovlll998l ). 

The discrepancy between the frequency scaling of 



core position observed in PKS 1257-326, and the typi- 
cal r cor e v~ x scaling observed in sources studied with 
VLBI may arise because different types of source lend 
themselves to different types of analysis. Only the largest 
core shifts can be detected with VLBI im a ging. In an 
analysis of 277 objects by iKovalev et al.l (|2008f h only 
10% gave reliable core shift measurements. The sources 
used in the detaile d, mult i- frequency VLB I study of the 
core shift effect by iSokolovskv et al.l ([20111 ) were selected 
based on their known large core shift, in addition to their 
bright optically thin jet features that were required to 
enable accurate registration of the images. 

More specifically, we suggest that the peculiar nature 
of the core-shift frequency dependence in PKS 1257-326 
is related to a number of other remarkable source prop- 
erties. The implied brightness temperature of the scin- 
tillating component in PKS 1257-326 is high, > 10 12 K 
(Bignall et al. 2006), and the source has exhibited such 
bright emission at least as long as IDV has been observed, 
since 1995. The persistence of this IDV over more than 
15 years is relatively rare amongst IDV sources. During 
the 4-epoch, one year duration MASIV survey, Lovell et 
al. (2008) found that sources which consistently exhibited 
IDV over all four epochs accounted for only 20% of all 
IDV sources observed. Moreover, the stability of the an- 
nual cycle in the source variability timescale (Bignall et 
al. 2003, 2006) suggests that the size of the source is re- 
markably stable on a timescale of years. Physically, this 
implies that the ultracompact jet in this source is both 
remarkably bright and stable. It is tempting to spec- 
ulate, therefore, that the hydrostatic confinement pro- 
vided by the strong pressure gradient suggested by our 
simple model is responsible for the observed stability of 
the jet. 

The Australia Telescope is funded by the Common- 
wealth of Australia for operation as a National Facility 
managed by CSIRO. The observations presented here 
were made by J AH as part of the CSIRO Astronomy 
& Space Science (CASS) Vacation Scholarship Program. 
JAH thanks Dominic Schnitzeler for assistance with the 
observing setup. Parts of this research were conducted 
by the Australian Research Council Centre of Excellence 
for All-sky Astrophysics (CAASTRO), through project 
number CE110001020. 
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APPENDIX 

DERIVING SOURCE STRUCTURE FROM THE TIME DELAY BETWEEN SCINTILLATION LIGHTCURVES 

Here we relate the time delay observed between the scintillations at two frequencies to the angular displacement 
between the centroids of the source emission at the two frequencies. The time delay depends not only on the velocity, 
but also on the axial ratio, R, and orientation of the interstellar scintillation pattern. 

For intensity fluctuations A/(r, t) measured at two locations r' and r' + r on the observer's plane at time t, the 
intensity correlation function is defined as, 

, (r , t)= (AW)A/ (; < + r.,- + ,)) (A1) 

We follow the treatment of Coles & Kaufman (1978), in which the simplest approximation is that the contours of equal 
spatial correlation comprise a family of similar ellipses, and the model for the spatial correlation function takes the 
form 

p(v,0) = f(^p], (A2) 
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where / is a monotonically decreasing function of |Cr| 2 . If the scintles are elliptical with axial ratio R and x axis is 
inclined at an angle j3 with respect to the major axis, then C takes the form 



C = 



cos 

— \/Rsmf3 \/RcosP 



(A3) 



The scintillation pattern is assumed frozen onto a screen that moves past the observer at velocity v, so there is a 
direction relation between the spatial and temporal dependence of the autcorrelation, 



Surfaces of constant correlation correspond to curves of constant \C(r — v£)| 2 . The maximum of p(r, t) occurs at a 
time lag At given by 

At = a • r (A5) 

where 

(C T C) • v 



(A6) 



Equations (|A2|) — ()A6|) are derived in Coles & Kaufman (1978), and are shown here for completeness. 

Now suppose a source possesses similar structure at two frequencies v\ and v<i, but that they are displaced by an 
angle 0. Then the scintilation fluctuations at v\, A/i(r, t), are identical to those at v<i, AJ2(r, t), except that they are 
displaced by a linear scale DO (see, e.g. Little & Hewish 1966): 

I 2 {v) = h{v - DO), (A7) 

where D is the scattering screen distance. The cross-correlation between l\ and I2 takes the form 

(AIi (I"', *') AJ 2 (r' +r,t' + t)> (A/j (r', t')Ah (r' + r- DO, f + 1)) 



(AhAI 2 ) {AhAI 2 ) 



p hh {v - D0,t). 



(A8) 



Onc e again, we assume that pi 1 i 2 is a montonically decreasing function of r that takes the form given by equation 
(|A2|) . Now, if the scintillations are measured at identical locations, we are interested in the max imum of pi 1 i 2 (—D0, t). 
In this case the time delay measured between I\ and I2 is the same as that given by eq. (|A5|) above, with r replaced 
by — DO, namely, 

At = -Da • = .^■v + (tf-l)(^xS)(vxS) a 

v 2 + (i? 2 -l)(vxS) 2 

where S = (cos /?, sin /3) is the direction along which the scintles are oriented. The cross product of a vector with S is 

the component of that vector that points orthogonal to the elongation axis. For instance, r x S is the component of r 
orthogonal to the long axis of the scintillation pattern. 
Finally, we note that the timescale of the scintillations can also be derived from the foregoing formalism, and is given 

by, 

&y/R 

4cint = , (A10) 

V + (i? 2 -l)(vxS) 2 



where cr, defined by eq. (|A2|) , is the overall scale factor of the scintillation pattern. 

CORE SHIFT IN THE PRESENCE OF A PRESSURE GRADIENT 

As discussed by lLobanovl (jl998l ). the frequency dependence of core posit ion may be influ enced by the pressure 
gradients in the external medium if hydrostatic confinement is important. lLobanovl ([ 19981 ) plotted the frequency 
dependence of core position as a function of the power-law index of the pressure profile, but neglected the effect of the 
changing Doppler factor along the jet. Accordingly, we present a derivation of the properties of the core (diameter, 
distance along the jet, and spectral index) for a simple model of a hydrostatically confined jet in the presence of a 
power- law pressure profile, accounting for the effect of a radially dependent Doppler factor. 

Let L be the path length through the source, d the diameter of the jet perpendicular to the jet axis, Da the 
angular size distance, the jet viewing angle, and a! v , the absorption coefficient in the source co-moving frame, at 
the rest frequency v' = ^-^v- Assuming the jet opening angle is small compared to the jet viewing angle (that is, we 



approximate the jet as a cylinder), we can approximate the optical depth as, 

d 



Slll(p 

cl v ,, (Bl) 



sin cj) 5 

since vol v is a relativistic invariant. For a power law electron distribution, the absorption coefficient is, 

l + z X"(«+ 2 - 5 ) 



cl v , = C(a)k e B a+1 - 5 {^-j^v 
where C(a) is a constant which depends only on the optically thin spectral index a. So now, 

A /I -I- r \-(«+ 1 - 5 ) 

Tv = ^-C(a)k e B a+1 - 5 ( i±£ ) . (B2 ) 

sm0 \ o J 

Hence, at the t v — 1 surface, which we identify with the scintillating component, or core: 



(a+2.5) 7 ^ ^ r r/^N d/^ 1 1-5+a 



i/y a+ ^ _ d(r) fe e (r) 



^0/ d(r )k e (r ) 



S(r)B(r) 



S(r )B(r ) 



(B3) 



ncxr 3np/4 



Tmin 0^ ^ 

npQ + 1.5) 



where ro is the position of the r v = 1 surface, or core, at frequency ^o- 

Hydrostatic confinement implies that the pressure inside the jet adjusts to the pressure of the surrounding medium. 
In that case, the diameter of the jet, d(r), is determined entirely by the pressure gradient of the external medium, 
and the lateral expansion of the jet dictates the run of magnetic field, particle density, and Lorentz factor along the 
jet. Consider a jet with ultra-relativistic equation of state, and relativistic particle energy distribution of the form 
N(j) = K e ~f~ a between some minimum and maximum Lorentz factor, 7 m i n and 7 max , confined by an ambient medium 
with pressure p ext oc r~ np . In that case, conservation equations imply that the following relations hold: 

d( r )(xr n *> /4 (B4) 

r(r)ar np/4 

£||Ocr- n ^ /2 
j B ±ocr -n P /4 

-3n p /4 
-n p /4 

k e ocr~ 

Following iKonigll (|1981[ ) and lLobanovl (|1998[ ), let us define k r such that 

r core ex (B5) 

and tub such that 

m B n p 

Boer — . (B6) 

In that case, ra# = 1 corresponds to a predominantly toroidal (perpendicular) magnetic field, while tub = 2 corresponds 
to a predominantly poloidal (parallel) m agne tic f ield . Approximating the radial dependence in Doppler factor as a 
power law of the form 5 oc r ns , Equations IB3I and IB4I give 

, K/4)[(1.5 + a)(2 + tub) - 1] - (1.5 + a)n s 

k r = — — . (B7) 

2.5 + a 

Note that n$ may be positive or negative depending on the initial Lorentz factor, T(ro), and jet viewing angle 0, and 
therefore the effect of the Doppler factor may be to steepen or flatten the relationship between r and v. 

For a self-absorbed synchrotron source, the flux density at the peak frequency (which we associate with the flux 
density at the observed frequency) is, 

5 m OC^ 2 ^ S -l/2^1/2. (Bg) 

This expr ession, when combined with the functions d(r), B(r), 5(r) and r(^), allows a prediction of the core spectral 
index (see Konigl 1981, equation 12). Here we define the spectral index of the core, a core , such that S v ex v acove , and 

5 1 r / ttib\ 1 
ow = 2 - [ n P y 1 + ~YJ + ns \ ' ( B9 ) 
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